Exact solutions and elementary excitations in the XXZ spin chain with unparallel 

boundary fields 
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By using a set of gauge transformations, the exact solutions of the XXZ spin chain with unparallel 
boundary magnetic fields are derived in the framework of the algebraic Bethe ansatz. In the easy- 
plane case, we show the elementary excitations are some kind of spinons without definite spin because 
the U(l) symmetry is broken, while in the easy-axis ferromagnetic case a spiral state is realized in 
the ground state. The correlation functions, the spin torque as well as the spin voltage for the later 
case are also derived. 



•4-* 



o 
o 



> ■ 
m ■ 

: 

(N ■ 

<n : 
o . 

-I— > ' 

a ■ 

a 

T3 ■ 

c : 
o . 
o ■ 



X 



The magnetic structure in a variety of condensed mat- 
ters has been an important issue in modern physics. In 
lower dimensions, the situation is especially distinctive 
and the spins may show a variety of exotic features due 
to the strong quantum fluctuations. For instance, the ele- 
mentary excitations in the one-dimensional (ID) spin liq- 
uids are doublet kinks M rather than triplet spin waves or 
magnons. Such quanta carry spin 1 /2 and are now known 
as spinons. Furthermore, in spintronics, spins may carry 
information so that extended capability and performance 
can be added into the electronic device. O Several ways 
to produce pure spin current (without electric current) 
have been proposed, such as spiral spin structure as a 
driving force, || spin Hall effect in ferromagnetic met- 
als, pj and ferromagnetic resonance. ||||(| However, so 
far no experimental realization of pure spin current has 
been achieved. The spiral structure is a long range or- 
dered state and therefore could hardly exist in ID mate- 
rials such as quantum wires. Some interesting questions 
one may ask are: Can two unparallel boundary magnetic 
fields generate a spiral spin state in ID? How does a spin 
liquid response to such fields? Since generally these fields 
break the U(l) symmetry of the system, what is the na- 
ture of elementary excitations in ID spin liquid under 
such fields? 

As proposed in Ref. [3] , the spin-supercurrent state may 
occur in easy-plane ferromagnets. In this work, we show 
that the spin structure, and consequently, the spin volt- 
age may be controlled by the magnetic anisotropy repre- 
sented by the boundary magnetic fields. By using a set of 
local transformations, the initial state for algebraic Bethe 
ansatz is established and the exact eigen-states of the 
model Hamiltonian are constructed. In the easy-plane 
case, the elementary excitations are almost the same as 
the usual spinons in spectrum but show spiral behavior in 
the real space and do not carry definite spin because the 
[/(l)-symmetry is broken. In the easy-axis ferromagnetic 
case, the ground state is a spiral state. The magnetiza- 
tion, the correlation functions, the spin torque as well as 



the spin voltage in the ground state are derived exactly. 
The model Hamiltonian we shall consider reads 

JV-l 

+ hi ■ &i + hi ■ <jn, (1) 

where a" (a = x,y,z) are the Pauli matrices, rj is the 
anisotropic parameter and hi = Ci(sin0,cos#, hi) and 
h-2 = ^(sinS, cos 6*, h^) are the boundary fields. If both 
h\_2 along the z-direction of the spin space, the system 
possesses U (1) symmetry, i.e., the z-component of the to- 
tal spin is conserved and the elementary excitations are 
the usual spinons. [Q. However, if one of the bound- 
ary fields has non-zero transverse component, the U(l) 
symmetry is broken and the spin structure may be dif- 
ferent. We note that the present model has a lot of im- 
portant applications in several physical systems. It has 
been providing valuable insights into non-equilibrium be- 
havior and complex dynamics |7|,p|,p|Jl0| and broken sym- 
metry. Interestingly, this model has been proven to 
be integrable p2j| via the open boundary transfer ma- 
trix method. |13| However, the spectrum has not been 
derived via the conventional Bethe ansatz scheme and 
leaves as an open problem in the field of integrable mod- 
els. The difficulty comes from the absence of an obvious 
initial state or pseudo vacuum (Without the boundary 
fields, it is the state where all spins polarized along the 
z-direction), thus the known Bethe ansatz can not be 
used directly. 

In the following text, by using a set of local gauge 
transformations, we shall construct an adequate initial 
state with which the eigen-states can be derived via the 
algebraic Bethe ansatz. To show our procedure clearly, 
let us review the integrability of the model (1) briefly. 
The integrability of the XXZ spin chain is associated with 
the Lax matrix 



1 



Lon{u) 



1 



2 sin 77 
1 

2 sin 77 
T o a n 



(sm(u + rj) + sinw) 
(sin(u + rf) — smu)<jQ<j 2 r 



(2) 



which satisfies the Yang-Baxter relation |14| , |15| 

R a pL an {u)Lf jn {v) — Lf3 n {v)L an (u)R al3 , (3) 

with R a p = L a p(u — v). The monodromy matrix Tq(u) = 
Lqi(u) ■ ■ ■ Lqn(u) also satisfies the Yang-Baxter relation. 
In addition, as demonstrated by Sklyanin [ fL3| , the so 
called reflection equation 

R a/ 3(u - v)K~(u)R a p(u + v)Kp{v) 
= Kp(v)R aP (u + v)K-(u)R al3 (u-v) (4) 

is also necessary for establishing the open boundary inte- 
grable models. Since the double-row monodromy matrix 



Uo(u)=T (u)Kc(u)T - 1 (-u) 



(5) 



also satisfies the reflection equation, the transfer matrices 
t(u) = troKQ (u)Uo(u) with different spectral parameters 
are mutually commutative, i.e., [t(u),t(v)] = as long as 
the dual reflection matrix Kq (u) satisfies the dual reflec- 
tion equation. []l3f Therefore, t(u) serves as the genera- 
tion function of the conserved quantities. For the XXZ 
spin chain, the general solution of K^(u) read |l^ , p^| 



Kq (u) = coo cos u — coi sin uuq 

+ sin(2u)(sin0<To +cos0cJq), 
(u) — -coo cos(u + ??) - coi sin(w + ij)<Tq (6) 
+ sin2(u + ri)(smd<TQ + cos OcTq), 

where coo, coi, coo, On, 9 and are arbitrary parameters. 
The Hamiltonian can be expressed as 



H = sin rj— ln£(u)| u= o — N cos rj, 
du 



(7) 



with C\ = — 2sin?7CQQ 1 , hi = — coi/2, C 2 = 2sin77c 00 1 and 
h 2 = -cm/2. 

An obvious initial state exists in the conventional al- 
gebraic Bethe ansatz scheme, acting on which, the lower 
left element of the 2x2 matrix Uo(u) is zero. However, 
in the present case, the initial state is not obvious be- 
cause Kq and Kq are generally not diagonal. To find 
a proper initial state, as used in solving the eight-vertex 
model, [ p^lJT^Jl^ l we introduce two sets of gauge matrices 
M n (u) and M n (u). Under the gauge transformations 

L 0n (u\m) = M~\ +m (u)L 0n (u)M n+m (u), 
S 0n {-u\m) = M~l m (-u)LQ^(-u)M n - 1+m (-u), (8) 

the lower left elements of Lo n and its inverse Son acting 
on a proper local state is zero. The above requirements 
define the forms of the gauge matrices as 



M n (u) = (X n (u),Y n (u)) . 
M n {u) = (Y n+ i(u), Y„_i(u)) . 



(9) 



where X n and Y n are 1x2 column matrices with upper 
elements e~ lu x n and e~~ m y n , respectively; their lower el- 
ements are defined as unity. The constants x n and y n 
read 

x n = -i e - ir i^+ n \ y n = -ie-"'^-"), (10) 

where [3 and 7 are two parameters determined by the 
boundary fields. In addition, we define 



X n (u) 



1 



-( e ,#n); 



•En Un 

Y n {u) = —^—(e m 1 -y n ), 
X n {u) = (-e ia ,x n ), 



Y n (u) 



X n - y n -2 
1 



-(e™,-y n ). 



(11) 



Xn+2 - Un 

In this case, the local vacuum of site n reads 

\u >™= x n+m \ t>„ +| |>„ . (12) 

It is easily to show that the lower left elements of 
*S'on(— u\m) and Lo n (u\m) acting on this state is zero, 
while the upper diagonal element of So n (— u\m) and 
the lower diagonal element of Lo n {u\m) shift \w >™ to 
\u >™~\ and the lower diagonal element of So n {— u\m) 
and the upper diagonal element of Lo n (u\m) shift \ui >™ 
to \lu >™ +1 , respectively. Now we can rewrite the trans- 
fer matrix as 



t(u) = tr Q K^(u\m )U (u\m ), 



(13) 



where 



K+(n\m\ - ( ^0 + ("Mll ^o + ( u Ml2 

~ { K+(u\m) 21 K+(u\m) 22 

Y m (-u)K+(u)X m (u) Y m (-u)K+(u)Y rn _ 2 (u) 
X m (-u)K+ (u)X m+2 (u) X m (-u)K+ (u)Y m (u) 



and 



Uo(u\m) 



A m (u) B m (u) 



C m {u) D m (u) 

Y„^ 2 {u)U{u)X m {-u) Y m (u)U(u)Y m (-u) 
X m (u)U(u)X m (-u) X m+2 (u)U(u)Y m (-u) 



(14) 



To ensure Kq (u|too) being diagonal, we parameterize the 
left boundary field as coo + zcoi = —2 cos(mo + 7)??, then 

t(u) = K^(u\m )nA mo (u) + K$ (u\m ) 22 D mo (u). (15) 

In addition, C m (u) must annihilate the global initial state 
\Q > m =n (S \uj >I™. This needs coo + ic i = ~2cos(9-d+ 



2 



r](j + N + m— 1)), which gives a constraint of coo, cqi, coo 
and coi. Under the above condition, |fi > m is a common 
eigen-state of A m (u) and D m (u), and B m (u) serve as the 
generating operators of the eigen-states of t{u) as usual. 
To derive the Bethe states, we need also the commutation 
relations between A m , D m and B n . With the help of the 
production relations of R\2{u\ — U2)X^ n+2 {ui)X^ n+1 (u 2 ) 
etc. [|l7| and the reflection equation of Uq(u), after some 
manipulation we obtain 

B m {ui)B m -2{u 2 ) = B m (u 2 )B m _ 2 (ui), (16) 
sin(ui + tt 2 ) sin(ui - u 2 - 7/) 
sm(ui - u 2 ) sm(ui + u 2 + rj) 
xB m (u 2 )A m (ui) 
sin(2?i 2 ) sin?7sin(ui — U2 — r\(m + 7+1)) 
sin(ui — u<x) sin 77(777 + 7 + 1) sin(277 2 + 7/) 
xB m (ui)A m (u 2 ) 
sin(— U\ — 112 + r](m + 7)) sin 77 sin 77 
sin 77(771 + 7 + I) sin(iti + U2 + 77) sin(2?j 2 + 7/) 

xB m (ui)D m (u 2 ), (17) 
D m+2 (ui)B m (u 2 ) = 
sin(ui — 7i 2 + 77) sin(tti + u 2 + 277) 
sin(ui - u 2 ) sin(iti + U 2 + 77) 

xB m (u 2 )D m (ui) 
sin(ui — U2 + 77(771 + 7+I)) sin(2ui + 277) sin 77 
sin 7/(777 + 7 + I) sin(7i! — tj 2 ) sin(2M 2 + 77) 

xB m (ui)D m (u 2 ) 
sin(2u2) sin(tti + u 2 + 77(777 + 7 + 2)) 
sin(iti + ?i 2 + 77) sin 77(777 + 7 + 1) 
sin(2iii + 277) 



sin(27i 2 + 77) 



B m ( Ul )A m (u 2 ), (18) 



where 



- sin(2zi + 77) sin 77(777 + 7 + 1) 
D m (u) = : ; -. ; — r D m (u) 



sin 77 sin 77(777 + 7) 
sin(?7 ( 777 + 7) + 2u + 77) 
sin 77(777 + 7) 



A m {u). (19) 



Assume the eigen-state of the transfer matrix take the 
form 

|$ >= B ma _ 2 {ui) ■ ■ ■ B mo _ 2M (u M )P > m , (20) 

with mo — 2M = 777. Acting the transfer matrix on this 
state and exchanging the position of A mo (u), D mo (u) and 
B mo - 2 {ui)B ma _ 4 (u 2 ) ■ ■ ■ B mo _ 2M (u M ) by the commuta- 
tion relations, we obtain two types of terms. One gives 
the eigen value of t(u) while in the other type one spec- 
tral parameter Uj exchanges with u. The latter terms are 
set to zero and we obtain the Bethe ansatz equations 



sinh 2JV (A 3 - \%n) 
sinh 2JV (A J + ±777) 



sinh(Aj + iai — \ir\) cosh(Aj + ia 2 — \ir\) 
sinh(Aj 
sinh(Aj 



iot\ + ^777) cosh(Aj — ia 2 + ^777) 
iai — 2 iil) cosh(Aj — 7a 2 — ^777) 



sinh(Aj + ia± + 5*77) cosh(Aj + 7a 2 + ^777) 



M 

n 

*+?■ 



sinh(Aj + A; — 777) sinh(Aj — A/ — 777) 
sinh(Aj + A/ + 777) sinh(Aj — A/ + 777) ' 



J 



1,2, 



,M. 



(21) 



The energy spectrum is given by 

sin 77 sin(ai — a±) 



E = (N - 1) cos 77 



smai sinai 



M 

£ 

3 



2 sin 77 



' cosh 2 Xj — cos 2 (^77) ' 



(22) 



with the paramctcrizations coo = —2 sin u\ cos a 2 , Coi = 
2i cosai sina 2 , cqo = — 2 sinai cos fi 2 , cqi = 



2? coscii sina 2 . Here we have put Uj 



77/2. The 



Bethe ansatz equations look very similar to that of the 
usual case. However, there is a restriction on the num- 
ber of "flipped" spins M. In fact, the two boundary 
fields are not independent with each other but related 
with a complex parameter 7 + mo and an integer M. 
That means given a right boundary field, the left bound- 
ary is restricted by the right one and M or vise versa. 
In addition, by exchanging the two fields, we can ob- 
tain another set of solutions with M — > N — M — 1 
and the same Bethe ansatz equations and spectrum. We 
choose coo + icai — — 2cos(t77o + n /)v and Coo + icoi = 
—2 003(6* — 6 + 7/(7770 + 7 — 1)) for real 77 (easy-plane case) 
and even N. In this case, M = N/2, N/2 - 1. This 
allows us to study the ground state and the low lying ex- 
citation properties simultaneously. In the charge neutral 
sector (M = N/2), the ground state is described by real 
Xj which fill the whole real axis in the continuum limit; 
the excitations are A — n-strings with 77-holes in the real 
axis. In the thermodynamic limit N — ► 00, the excitation 
energy reads [fill 



e(A^--A£)=£>(A?)=£ 



7r sm?7 



J_ 

27) 



77 cosh 



(23) 



i.e., the spectrum takes the same form of the case with- 
out the boundary fields, where A^ is the position of j-th 
hole in the spectral space and e(A^) is the energy of a 
single spinon. |l{| The other type of excitations can be 
constructed by considering the process of M = N/2 to 
M = N/2 — l, i.e., a two- hole state based on the ground 
state. The excitation energy simply reads 



e(X h 1 ,X h 2 ) = e(X h i) + e(X%), 



(24) 



where Xi 2 are the positions of the two holes in the spec- 
tral space. This state corresponds to the spin triplet exci- 
tation in the usual case. The above discussion shows that 



3 



the elementary excitations are not changed by the bound- 
ary fields in spectrum, i.e, they are also described by 
spinons formally (additive spectrum). Since the spinons 
are some kind of topological excitations, it is rather 
strange that the broken symmetry does not alter the 
spectrum. However, the spin structure in the real space is 
indeed changed by the boundary fields. The initial state 
is in fact a spiral state with the spins rotate site by site 
in the x — y plane. The "spin- wave" generators B m (u) do 
not "flip" the spins on different sites to an uniform direc- 
tion. Therefore the spins have winding behavior around 
the z-direction. 

For the easy-axis ferromagnetic case, we put rj — ir + 
i\ with x rea l an d study the coo + icoi = — 2cos((7r + 
ix){ m o+l)), and c o + «coi = —2cos(6—6+(ir+ix){mo+ 
7 + N — 1)) case. The ground state is described by the 
reference state |fi > m °. Several quantities in the ground 
state can be derived easily. For instance, the z-direction 
magnetization is 



<M z >=^ + 0(N- 1 ), 



(25) 



implying that the boundary fields change the magnetiza- 
tion only in order of l/N. The correlation functions can 
also be derived as 



< o z n a z m >= tanh/(n) tanh/(m) 



(-!)«+"' cos 2 (g) 



<0-„°W»>- cosh cosh /( 7 



(26) 



T VrrV 



(-1)™+™ siiU( 3 ) 
cosh /(n) cosh /(m) 



with f(n) — (Re(mo + 7) + n — l)x + 7r/m(mo + 7) and 
g = nRe(mo + 7) — x/m(rao + 7) — Obviously, the z-z 
correlation function tends to 1 when both m and n tend 
to infinity, indicating the long-range order of the ground 
state. Define the quantum spin torque operator as S n — 
jff n x <7„+i- The expectation value of this quantity in 
the ground state reads 



< si > = (-iy 

< s* > = (-iy 

< s* > =0. 



, sm(g) sinh[/(n) + y/2] cosh(x/2) 

2 cosh f(n+ 1) cosh /(n) 
..+1 cos(g) sinh[/(n) + y/2] cosh(x/2) 
2 cosh/(n + 1) cosh /(n) 



JV-l 



The spin voltage can also be calculated as V" = Yl' n ={ < 
S% >, a — x, y, z. 

In conclusion, we develop a method to derive the eigen- 
states of the XXZ spin chain with unparallcl boundary 
fields. For some special choices of the boundary fields, 
the ground state and the low lying excitations for the 
easy-plane case are obtained exactly. The spectrum of 
the excitations is of the same form as that of the usual 
spinons without definite charge (spin), suggesting that 



there exist some hidden symmetries in the system, an 
interesting question needs to be further explored. For the 
easy-axis ferromagnetic case, the ground state is a pure 
spiral state which possesses non-zero spin torque in the 
x-y plane. This provides a possible scheme to generate 
driving force for spin current in quasi ID systems. 
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